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Abstract. 


We  introduce  a  new  class  of  multivariate  new  better  than  used  (MNBU) 
life  distributions  based  on  a  shock  model  similar  to  that  yielding  the  Marshall- 
Olkin  multivariate  exponential  distribution.  Let  T^,  . ..,  be  indepen¬ 
dent  new  better  than  used  (NBU)  life  lengths.  Let  F(t^,  ...»  t  )  be  the 
joint  survival  function  of  min  T.,  i  “  1 ,  ...»  n,  where  A^,  ....  A^  are 

j*A  J 

n 

nonempty  subsets  of  {1 . M}  and  n  A^  =  {1,  ...,  M)  .  F(tj,  •••»  tn )  is  said 

i*  1 

to  be  a  MIIBU  survival  function.  Basic  properties  of  MNBU  survival  functions 
are  derived.  Comparisons  and  relationships  of  this  new  class  of  MNBU 
survival  functions  are  developed  with  earlier  classes. 


1.  Introduction 


^  The  univariate  *new  better  than  used*  (N3U)  class  of  life  distributions 

was  sliown  by  Marsliall  and  Proschan  (  1972  )  to  play  a  key  role  in  the  study 

of  maintenance  policies.  See  also  Harlow  and  Proschan  (  1975),  Chap.  6. 

,T 

.■>In  this  paper,  -we  introducers  multivariate  version  of  the  NBU  distribution 
based  on  a  physical  model.  Shocks  occur  in  time  which  cause  the  simultane¬ 
ous  failure  of  subsets  of  components.  The  interval  of  tine  until  the  occur¬ 
rence  of  a  shock  destroying  a  given  subset  of  components  is  governed  by 
an  NBU  distribution.  The  occurrence  times  are  mutually  independent.  Note 
tlie  similarity  of  this  model  with  the  shock  model  leading  to  the  tlarshall- 

Olkin  ( 1967 )  multivariate  exponential  (MVE  )  distribution.  In  the  Marshall- 

' i  <7 

Olkin  model  shock  times  have  exponential  distributions;  in  .our  model,  shock 
times  have  NBU  distributions. 

rv. 

Other  versions  of  the  multivariate  NBU  distributions  have  been  intro¬ 
duced  and  studied.  See,  e.g.  ,  Marshall  and  Shaked  (1979a,  b ).  Our  model 
may  be  of  interest  in  certain  applications,  since  the  underlying  notion 
derives  from  a  shock  model. 

In  Section  2  we  give  two  equivalent  formulations  of  the  MNBU  class 
and  obtain  its  properties.  In  Section  3,  wa  consider  other  classes  of  multi¬ 
variate  new  better  than  used  life  distributions  and  compare  them  with  our 


MNBU  class 
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!.  Definitions  and  Properties. 

In  this  section  the  : IIIBU  class  is  defined  and  its  properties  are  studied: 
A  fatal  shock  model  generating  distributions  in  this  class  is  formulated. 

This  model  is  a  direct  generalisation  of  the  Marshall  and  Olkin  (1976)  MVE 
model . 

We  begin  by  giving  two  equivalent  definitions  of  the  IINBU  distributions: 

Definition  2. 1.  A  random  vector  T  =  (T  ,  . ..,  T  )  is  said  to  be  a  HMBU 

random  vector  if  it  has  a  representation  T.  =  min  T,,  where  TA.  Ac  I  , 

“  i  .  A  A 

ieA 

are  independent  NBU  rand  on  variables  (possibly  dengerate  at  0  or  00 )  and  I 
is  the  class  of  nonempty  subsets  of  (1,  n). 

Definition  2.2.  A  random  vector  _T  5  (Tj,  ...,  T  )  is  said  to  be  a  IINBU  random 

vector  if  it  has  a  representation  T  =  min  X.,  where  X  ,  ...,  X  are  inde- 

1  .  c  J  1  m 

16  i 

pendent  NBU  random  varibles  (possibly  degenerate  at  0  or  <*)  and  4  ^  c 

n 

{ 1 ,  •  •  •  ,  11}  ,  i  ~  1 ,  a ,  a nd  u  B { 1 ,  » .  •  ,  m). 

1=1 

Remark  2. 3.  The  two  equivalent  formulations  above  permit  us  to  use  whichever 
is  more  convenient. 

In  Definition  2.1,  when  the  random  variables  T^,  A  e  I  ,  are  exponential, 

then  jT  is  the  MVE  rcndon  vector.  Recall  that  _T  can  be  viewed  as  the  vector 

of  life  lengths  of  n  components  subject  to  fatal  shocks  from  independent 

sources.  For  every  A  e  J,  T  is  the  random  time  at  which  a  shock  occurs 

A 

which  simultaneously  destroys  all  the  components  whose  indices  form  the 
set  A. 

Let  F(t,,  ...,  t  )  *  P(T,  >t,,  . . . ,  T  >t  )  be  the  joint  survival 
l  n  11  n  n 

function  of  T^,  where  T  is  MNBU.  Eq.  (2.1)  expresses  F(tj,  ...»  t^) 


\ 


in  terms  of  F.  A  e  I,  where  F  is  the  survival  function  of  T  : 

*  A  A 

(2.1)  F(  t . .  )=  n  F (max  t .  ),  t.  ^  0 ,  i  =  1,  n. 

n  AcI  A  ieA  1  1 

The  following  lemma  shows  that  F(tj,  t  )  enjoys  a  property  similar 

to  the  defining  property  of  UB U  random  variables. 

Lemma  2.4.  Let  F(t  ,  t  )  be  defined  bv  (2.1).  Then 

- 1  n 

(2.2)  F(tj  +  s,  ...,  ta  +  s  )  <  F(tlS  ...,  tn)F(s . s)  for  all 

s  ^  0,  t ^  ^  0,  i  =  1 . . 

Proof .  Since  max(t  +  s)  =  max  t.  +  s,  and  T  is  NBU  for  each  A  e  I ,  we  have: 
ieA  1  ieA  1 

F(t  +s,  ...,t  +s)=  II  F  (max(  t  +  s  ) )  <  II  F  (max  t  )F  (s  )  ■ 

1  n  A el  A  ieA  1  Ael  A  ieA  1  A 

F(  t  ,  « « «  ,  t  )F(s,  . .  .  ,  S  ).  j  I 
1  n 

Remark  2. 5.  Note  that  (2.2)  can  'ye  expressed  as  P(T^  >  t}  +  s,  .  .., 

T  >  t  +  s|T.  >  s,  ...,  T  >  s)  <  P(T,  >  t.,  ....  T  >  t  ).  This  asserts 
that  the  joint  survival  probability  of  n  components  each  of  age  s  is  less 
than  or  equal  to  the  joint  survival  probability  of  n  new  components.  Another 
alternative  interpretation  of  (2.2)  may  be  obtained  by  rewriting  it  as 

P(T1  >  C1  +  s’  *•*  *  Tn  >  Cn  +  S,T1  >  V  '  Tn  >  tn)  -  P(T1  >  S . Tn  >  S  ) 

This  states  that  a  series  system  of  n  components  of  ages  t^»  . ..,  t^  is 

stochastically  shorter-lived  than  is  a  series  system  of  n  new  components. 

Remark  2.6.  A  multivariate  new  worse  than  used  (MNWU )  random  vector  T  can 

be  defined  as  in  Definition  2.1  (Definition  2.2)  where  now  T^,  Ael, 

(X^,  i  ■  1,  ...,  M)  are  assumed  to  be  independent  11WU  random  variables.  If 

F(t.,  . ..,  t  )  denotes  the  joint  survival  function  of  T  ,  . ..,  T  ,  then 
In  in 

we  can  easily  show  that  F(t.  +  s.,  ....  t  +  s  )  >  F(t,,  ...,  t  )F(s.,  s  ) 

li  n  n  —  i  ni  n 
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Note  that  in  the  Lt.’UU  case,  the  s  values  may  differ,  while  in  the  MNBU  case, 
the  s  values  must  be  the  sane. 

The  following  lenma  establishes  bounds  for  the  joint  distribution  and 

the  joint  survival  function  of  MN3U  random  vectors. 

Lemma  2.7.  Let  T  =  (T, . T  )  be  MNBU  and  let  F(  t t  )  and 

-  —  l  r.  In 

F(t . ,  ...,  t  )  be  the  ioint  distribution  and  tiie  joint  survival  function 
In 

of  T  ,  ...,  T  respectively.  Then 
1  n  n 

(i)  K(tt,  ...,  t  )  >  ii  (  ..  F  ft  ;), 

ic A 
n 

(ii )  F(  t  ,  ... ,  t  )  >  nil  -  n  F  (t.  ;]. 

1  n  "  i-i  i«A  A  1 

Proof .  Since  ,  ...,  T  are  increasing  functions  of  independent  random 
variables,  they  are  associated.  The  results  in  (i)  and  (ii)  follow  readily 
from  well  known  inequiilies  for  associated  tando.n  variables.  II 

The  following  theorem  shows  ‘•.hat  the  Ni.'HU  class  lias  nany  desirable 
propc  rt ies. 

Theorem  2.8.  The  following  properties  hold  for  the  MHBU  class: 

(PI  )  Let  T  be  an  NBU  randon  variable.  Then  T  is  1-dimensional  MNBU. 

(P2)  Let  Tt  ,  . . ,  T  be  independent  N8U  random  variables.  Then  T^  is  MNBU. 

(P3 )  Let  T  he  MNBU.  Then  (T.  ,  ...,  T  )  is  k-dinensional  MNBU , 

™  1  «  X . 

I  K 

1  ^  i.«  K  •••  K.  i.  ^  n .  k  “  1 ,  •  o  •  •  n  * 

—  a  k  — 

(P4  )  Let  T  be  MNBU  and  T*  *•  min  T  ,  i>  jt  B  c  { 1 ,  . . .  ,  n)  ,  j  -  1,  ...,  m. 

d  icd .  d 

J 

Then  T*  is  MNBU. 

(P5)  Let  T  be  MNBU  and  a  >0,  1*1,  ...,n.  Then  min  a.  T  is  NBU. 

l<i<n  1  x 

(P6)  Let  T_  be  n-dimeroional  MNBU,  T'  be  m-dimensional  MNBU,  and  T,  T*  be 
independent.  Then  (T,  T '  )  is  (n  +  n  )  -  dimensional  MNBU. 

(P7  )  Let  T  be  tUhlU  and  let  t  be  the  life  function  of  a  coherent  system. 


Then  t(T)  is  NBU 
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'P3)  Let  3:  (0,  «)  *  [C,  ~)  be  a  strictly  increasing  function  such  that 

g(x  +  y)£gfx)  +  g(y  )  for  all  x,  y.  Let  T[  be  MNBU ,  then  T '  i  (g(T^  ) . g(T^  )  ) 

is  MNBU. 

Proof .  (PI)  and  (P2)  are  obvious. 

(P3)  and  (P4):  Since  ( P3  )  is  a  special  case  of  ( P4  )  we  need  only  prove  ( P4  ). 


Let  T.  -  r.iin  X  ,  i  =  1,  n.  Then  T*  =  min  X  where 


&S  ' 


S'.  =  u  S  ,  j  =  l,  n,  and  thus  T*  is  MNBU. 

3  iCB.  1 
J 

(P5)  Let  T.  =  min  T  ,  i  -  1  ....  n.  Then  min  a.T.  «  inin{(min  A  )T  }  ,  an 

i£  A  ni<_i£n  A£  I  1 f  A 

NBU  random  variable,  since  a  series  system  f  independent  :!BU  random  vari¬ 
ables  is  NBU. 

(PC)  The  proof  is  obvious. 

fP7)  Let  t(T)  ■  max  min  T  ,  where  P  ,  ...,  P  are  nonempty  subsets  of 
l<r<P  i€Pr  1  1  r 

{ 1 ,  . . .  ,  n}  .  But  T  =  nin  X.,  6  t  S .  c  (1,  ...,  11)  ,  i  =  1,  ...,  n.  Thus 

j«si  3 

r (T_)  H  r  '(>. )  =  max  min  X  ,  where  A  **  U  S  ,  r  *  1 ,  ...,p.  Since  a 

l<r<p  ieA  r  i€P 

-  r  r 

coherent  system  of  independent  NBU  components  has  NBU  life  length,  the  desired 
result  follows. 


(P8)  Let  T  **  min  X.,  ^  ji  s  c  (l^  ...,  fl) .  Since  g  is  increasing,  we 

J«si  3 

have  g(T. )  ■  min  g(X, ),  i  =  1,  ...,  n.  Clearly  g(X. ),  ...,  g(X  )  are  inde- 
jeS1  3  " 

pendent  NBU  random  variables  and  consequently  g(T^),  ...,  gd^)  is  MNBU.  II 

l/e  conclude  this  section  by  giving  various  necessary  and  sufficient 

conditions  for  an  MNBU  random  vector  to  be  MVE. 
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Theorem  2.9.  Let  be  HNBU.  Then  the  following  conditions  are  equivalent: 

(i  )  T  is  MVE. 

(ii)  rain  a  T  is  exponential  for  all  a  >  0,  i  =  1 . . 

l<i<n  1  1 

(iii)  T^  has  exponential  minimums. 

(iv)  T  is  exponential  for  i  =  1,  n. 

(v)  niin  T  is  exponential. 
l<i<n  1 

Proof.  It  suffices  to  show  that(iv)  — >  (ii)  and  (v)—“ ■>  (i).  We  only  prove 

that  (v)—>  (i)  since  the  proof  of  (iv)*“— >  (ii)  is  similar.  Let  T  *  rain  X., 

jcSi  J 

b  t  c  {1,  •••,  II ) ,  where  Xj ,  ...»  X,;  are  independent  NBU  random  variables. 

flow  min  T  *  nin  X  ,  which  is  exponential.  Consequently,  each 
l<i<n  l<j<m  J 

X^  is  exponential,  and  so  is  MVE.  II 

3.  Other  Classes  of  Hultvariate  New  Better  than  Used  Distributions  and  Their 
Relation  to  the  tlNBU  Class. 

Several  alternative  definitions  are  available  of  multivariate  life 
distributions  extending  the  univariate  concept  of  I1BU.  Each  of  these  classes 
satisfies  some  of  the  properties  which  one  would  expect  for  a  class  of  mul¬ 
tivariate  new  better  than  used  distributions.  In  this  section  we  compare 
the  MNVW! class  with  some  of  these  other  classes. 

Consider  nonnegative  random  variables  T, ,  ...,  T  whose  joint  distribution 

i  n 

satisfies  on  of  the  following  conditions: 

(A)  T^,  ...,  are  independent  and  each  is  an  NBU  random  variable. 

(B)  (T. ,  . . . ,  T  )  is  MNBU. 

l  n 

(C)  For  all  a.  >  0,  i  *  I,  ...,  n,  min  a.T  is  NBU. 

l<i<n  1  1 


( [) )  For  each  d  /  A  c  {1,  ...,  n},  nin  T,  is  ai  ‘  Ti  l'  random  variable. 

ieA  1 

(E)  Each  T^  is  an  IJIJU  random  variable. 

Each  of  the  classes  of  multivariate  distributions  defined  by  (A)  -  (E) 
may  be  designated  as  a  class  of  multivariate  new  better  than  used  distri¬ 
butions.  h'e  now  compare  these  classes.  Clearly  (A)  — >  ( 'i )  —■>  (C)  — >  (D)  — >  (E 
The  following  examples  (see  Esary  and  marshal  j.,  1974  )  show  that  no  other 
implication  among  the  above  classes  is  possible. 

Example  3. 1 .  Let  T,  =  nin(l!,  57),  T0  =  nin(V,  '.'),  where  U,  V,  i'  are  inde¬ 
pendent  exponential  random  varibles  with  parameters  Xj  =  ,'  =  =  1. 

Then  (T^,  T^)  is  MI7BU,  but  T^  ,  T^  are  not  independent.  Thus  (A). 

Example  3.2.  Let  T  j  -  2Tj ,  T^  =  ,  where  Tj ,  T^  are  defined  in  Example 

3.1.  Obviously  min(a^T|,  a^T^)  is  fTBIJ  for  all  a^  >  0,  a^  >  0.  However 
(Tj,  T7  )  is  not  iiVE.  By  Theorem  2.7,  (Tj,  Tp  is  not  I1NBU.  Thus  (C)  -+>  (B). 
Example  3. 3.  Let  T ^  ,  T^  be  as  in  Exanple  3.1  and  let  (T*,  T^ )  =  (min(U,  57),  \\1). 
Let  F(tj,  t2)  =  pFT  T  (tj,  t2)  +  (1  “  p)FT*  T*(  1 1  >  t7),  where  0  <  p  <  1. 

Let  (Tj,  Tjj)  be  the  bivariate  random  vector  whose  joint  survival  function 
is  F(tj,  t2)>  Obviously  Tj,  Tjj,  and  min  (Tj,  T,j )  are  exponential,  but 
min(liTj,  T^  )  is  not  f!BU.  To  see  this,  let  F(  t  )  =  P(nin(5iT  j,T’  )  >  t)  = 
pe  +  ( 1  -  p  )e  ^,t.  It  is  easy  to  verify  that  F(2t)  >  [  F(  t )]  for 
sufficiently  large  t.  Thus  (0)“jt>  (C). 

Example  3.4.  Let  l),  V,  and  57  be  as  in  Exanple  3.1.  Leu  F(tjf  t2)  “ 

pFu  v(tl*  C2^+  ^  ”  P)FVJ  0^C1,C2^’  where  0  <  P  <  It  is  easy  to  verify  that 
F(tj,  t2)  is  the  joint  survival  function  of  a  bivariate  random  vector  whose 
marginals  are  NBU  but  whose  minimum  is  not  IJBU.  Thus  (E)^>  (D). 


Esary  and  .'larshall  (1974)  show  that  if  T.  has  exponential  ninimums, 
then  there  exists  an  liVE  random  vector  T*  such  that  t(T)  and  t(T_*  )  are  iden¬ 
tically  distributed  for  all  coherent  life  functions  x.  Unfortunately  the 
class  (D)  above  does  not  enjoy  this  property,  as  is  illustrated  by  the  follow 
in£  example. 

Example  3. 5.  Let  T^,  'l  ^  be  independent  exponential  random  variables  with 
parameters  A^  >  0,  >  0,  respectively.  ThenCT^vT^i  )  has  IffiU  minimums. 

Now  assume  there  exists  (Tj,  T^)  =  (nin(U,  U),  r,in(V,  W)  such  that  t(T) 
and  t(T' )  are  identically  distributed,  where  U,  V,  and  W  are  independent 
IIBU  random  variables.  We  then  have  F  •  G  °  H  =  and  F  •  H  =  F.,  where 
Fj,  F2,  F,  G,  and  11  are  the  survival  functions  of  TjVr^,  T.,,  U,  V,  and  U 
respectively.  This  leads  to  the  conclusion  that  C  =  F^/F^,  which  is 
impossible. 

Finally,  we  present  two  additional  classes  of  multivariate  new  better 
than  used  distributions  and  compare  them  with  the  MNBU  class.  The  first 
of  these  two  classes  is  due  to  ’larshall  and  Shaked  (  1979a),  the  second  is 


essentially  due  to  Block  and  Savits  (1979). 

Definition  (F).  A  random  vector  T  is  said  to  be  multivariate  new  better 
than  used  (F)  if  P^  >  0 ,  i  -  1 ,  . . . ,  n}  -  1  and  P(T  e  (a  +  B  )A>  < 

P{T  e  aA}P{T  e  tfA}  for  every  a  >_  0 ,  S  >  0,  and  every  open  upper  set  A  c  [0,  ")n. 


Definition  (G).  A  random  T  is  said  to  be  multivariate  new  better  than  used  (G) 
if  T^  has  a  representation  T  =  },  X  ,  where  X  ,  ...,  X  are  independent 

jcSj  J 

NVU  and  6  +  c  {1,  . ..,  Il),  i  ■  1,  ...,n. 

The  following  lemma  shows  that  the  MNBU  class  contained  in  MNBU  (F ). 
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Ler.ina  3.6.  Let  T  be  MMBU.  Then  T.  satisfies  the  conditions  of  Definition  (F). 

Proof.  Obviously  P(T.  >0,  i  =  1,  . n)  =  1.  Mow  T,  *  rain  X.,  where 

1  jeSj  1 

X, . X  are  independent  TIB U  randon  variables,  and  t>  t  c  {1,  . , 

i  =  1,  n.  The  desired  result  follows  immediately  by  Property  3.4  of 

Narstiall  and  Shaked  (  1979a).  !| 

The  following  examples  stow  that  no  other  implication  holds  between 
our  MMBU  class  and  the  MMBU  (F )  or  the  HNBU  (G)  classes. 

2 

—  "**?  ^  -^v 

Example  3.7.  Let  F(x,  y)  =  e  }  x  0,  y  0.  It  can  be  shown  that 
the  bivariate  random  vector  (X,  Y)  whose  joint  survival  function  is  F(x,  y) 
satisfies  (F).  Theorem  2.7  sliows  that  (X,  Y )  cannot  be  T1TIBU. 

Example  3.B.  Let  U,  V,  and  U  be  independent  exponential  random  variables 
with  parameters  Xj  ?  \„  and  X ^ 2  >  0  respectively.  Let  T}  =  min(U,  U)  and 
T2  -  nin(  V,  U).  Clearly  (T^ ,  T2>  is  IIMBU,  but  it  is  not  MMBU  (C).  To  see 

this,  assume  T^  =  X  +  Z ,  T2  =  Y  +  Z ,  where  X,  Y,  and  Z  are  independent  IJBU 

random  variables.  Since  T^  is  exponential,  it  follows  that  either  X  is 
exponential  and  Z  degenerate  at  0,  or  vice  versa;  similarly  for  Y  and  Z. 
Consequently,  T^  and  T2  are  either  independent  or  identically  distributed, 
which  is  impossible. 

Example  3.9.  Let  X,  Y,  and  Z  be  independent  with  absolutely  continuous 
distributions.  Let  T  -  X  +  Z,  T  =  Y  +  Z.  Then  (Tj,  T2  )  is  MMBU  (G), 

but  cannot  be  IIMBU.  For  if  (Tj,  T2 )  were  MMBU,  then  T^  and  T2  would  be 

independent,  which  is  not  the  case. 

Remark  3. 10.  In  Example  3.9,  observe  that  (Tj,  T2  )  ■  (X,  Y)+  (Z,  Z). 

This  shows  that  the  MNBU  class  is  not  closed  under  convolution. 
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